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Abstract, l^e construct almost invariant subspaces and the corresponding effective 

Hamiltonian for magnetic Block hands. We also discuss the question of the dynam- 

^ ics related to the effective Hamiltonian. We assume that the magnetic and electric 

Q ! potentials are slowly varying perturbations of the potential of a constant magnetic 

O ' field and a periodic lattice potential, respectively. 
CN ! 

^ ■ 1. Introduction 

In [5] we constructed wave packets for adiabatic perturbations of Schrodinger 
operators in periodic media. The recent work of Panati-Spohn-Teufel, [19], led us 
to consider the relation of those constructions to effective Hamiltonians. In §3 of 
this article we give a simple derivation of effective Hamiltonians for these problems. 
■ The main simplification in our method is the omission of the Floquet-Bloch 

transformation. This transformation has many nice properties. In particular, it is 
unitary, and this makes in useful in studying spectral properties of operators. In the 
O \ work of Helffer-Sjostrand [9] and Gerard-Martinez-Sjostrand [7] this transformation 
^ \ was used quite effectively in the computation of spectra, both of perturbed and ef- 
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fective Hamiltonians. However, if one is simply interested in effective Hamiltonians, 
Ph' the Floquet-Bloch transformation requires that one transform the Hamiltonian by 
r-| ■ a Fourier integral unitary operator only to transform it back at the end of the 
\ calculation. In this article we need to assume that eigenspaces of the unperturbed 
^ \ Hamiltonian depend smoothly on quasi-momentum, and form trivial bundles over 

a fundamental domain for the dual lattice. 

There is also the interesting question of how one interprets the lower order terms 
in the effective Hamiltonian. If one considers the propagation of observables in the 
5^ \ Heisenberg picture, it is natural to think of these terms as lower order corrections 
to the dynamics. This point of view is adopted in [19], and it is implicit in [3] and 
[4]. However, the highest order contributions of these terms to the wave packets are 
in a precession of the phase. Thus in [5] we did not include them in the dynamics, 
and did not see how to reconcile our results with those of [3] and [4] . It now appears 
that the two points of view complement each other instead of conflicting. 

2. Preliminaries 

The Hamiltonian for an electron in a crystal lattice V in M"^ in the presence of a 
constant magnetic field u = {uji^uj2tOJ^) is given by 

1 I d uj X X 



where y is a smooth, real-valued potential, periodic with respect to F. Here m and 
e are the mass and charge of the electron. To simplify notation we will use units in 
which h = 2m = e = 1. 

We will assume that F is generated by the basis {ei, 62, 63} for R^, 
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and let E be the fundamental domain {Y^^=i ij^j^ ^ [0; !)}■ We will use the dual 
lattice r* = e^Z + + 632, where e| • = 27r5jk, with the fundamental domain 

^* = {Ef=i*ie*,t,e[0,l)}. 

To realize Hq as a self-adjoint operator in L^(M^) we define it first on the 
Schwartz functions iS(M"^), and then take the Friedrichs extension. The resulting 
operator commutes with the magnetic translations introduced by Zak [24], 

T^/(a;) = e^<->^-'^>/V(^-7) (3) 

for 7 e r. We assume that 

(w, r X r) C 47rZ. 

With this assumption G = {T^,7 G F} is an abelian group, and we can reduce Hq 
by the eigenspaces of G, i.e. setting 

Vk = {ue HlM^), T^u = e-'^-^u, 7 e F}, (4) 

considered as a subspace of LP'{E), Hq restricted to is self-adjoint with compact 
resolvant. We denote its spectrum by 

E^{k)<E2{k) < ... 

Then by standard results the spectrum of Hq as an operator in L^(M^) is equal to 

UfeeE* VJ'^^^EUk). 

Note that, since Vk+'y* = Vk for 7* G F*, Emik + 7*) = Em{k). 

Standard perturbation theory shows that the function Em{k) is continuous for 
/c G M"^ and real analytic in a neighborhood of any k such that 

Em-i{k) < Em{k) < Em+i{k) (5) 

The closed interval = UkeE*Em{k) is known as the "m-th magnetic Bloch 
band" in the spectrum of Hq. 

In what follows it will be convenient to replace Hq acting on Vk by 

HQ{k) = e-'^--HQe'^-- = f-i-^ + + /cV + V{x) 



V = {ue Hl^{R^), T^u = 7 G F}. 



with the domain 

V = {ue i?iocV- />-7^ 

for all k. As with I>fe, we consider P as a subspace of L^{E). 

Assumption A. For a given rn we will assume that E^ satisfies (5) for all k. 

Under this assumption we can choose the eigenfunction k) associated to 
Emik) to be a real-analytic function of k with values in D, such that 

Hoik)^ ik) ^ Em{k)^{k) lor &\\\ [ \^(x,k)\^dx = 1. 
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Assumption B. We assume that 

*(a;, /c + 7*) = e^T*-^*(a;, k), 7* e T*. 

This assumption makes the complex hue bundle of the eigenspaces a trivial 
bundle over the torus , M'^/F*. In general one has 



*(a;, /c + 7*) = e*(^*-^+^^'='^*»*(a;, k), 7* e T*, 

where 9{k, 7*) is real-valued, and determines the structure of the eigenspace bundle. 
Since 

6{k, miel + 771262 + 777-363) = mid{k, el) + m20{k, 62) + 'msO{k, 63) mod 27r, 

when 9{k,j*) is nonzero, the derivatives of ^' with k will be unbounded and 
will not belong to the class of symbols B which we introduce below. Thus we need 
Assumption B. 

Remark 1. The general method of constructing effective Hamiltonians which 
we give here will apply under the weaker hypothesis: for a given m there exist p and 
q such that 

Em-p-i{k) < Em-p{k) and Em+q{k) < Em+q+i{k) for all k. 

However, in this case the effective Hamiltonian will he a matrix operator acting on 
functions with values in C^+^+-'^, as in [6], [7], [9] and [19]. 

3. Main Result 

The adiabatically perturbed Hamiltonian is 

o \ 2 

O iO X X 



= \^-i— + + A{ex) j +V{x) + W{ex), 

where W and A — [Ai, A2, A^) are smooth, and bounded together with all of their 
derivatives. As before, we define first on >S(M^), and then take the Priedrichs 
extension to get a self-adjoint operator in L^(]R^). 

The essential step in applying multi-scale techniques is simply to consider y = ex 
as a new independent variable in H^. Let 

Then, for u{x, y) we can define w{x) = u{x, ex) and conclude that 

[H^u]{x,ex) ^[H^w]{x). (6) 

The identity (6) enables us to solve the Schrodinger equation for uniformly in 
e by solving the Schrodinger equation for if^ uniformly in (y, e). The latter might 
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Let B denote the subspace of C°°(M'^ x R'^ x R-^) consisting of functions of the 
form 

P{x, y, k, e) = Po{x, y, k) + ePi {x,y,k)-\ h Pn{x, y, k) 

such that P{x + 7, y, k, e) = e*<'^^^'^>/2p(a;, y, k, e) and 

sup \\d^d^Pj{;y,k)\\L^^E) < 00, for all a,pe N\ 
y,k 

To P e B we associate the e-pseudo-differential operator 

P{x, y, eDy, e)f{x, y, e) = (27re)-3 J e^'= (^--)/^P(x, y, k, e)f{z)dzdk, f e <S(R^). 

Note that here we are using the standard quantization - as opposed to the Weyl 
quantization. Our main result is the following: 

Theorem. For every N E N there exist P^ = Fq + eFi + ■ ■ ■ + F^ € B and 
^eff — ^0 + ^hi + • • • + e^/ijv G i3 (independent on x) such that 

H,{Pn{x, y, eDy, e)u) - Pn{x, y, tDy, t)H^ff{y, tDy)u = 0(e^+i) (7) 

for u e <S(M^). Moreover, considered as an operator from L^(R^) into L'^{E x R^), 
Pjv is approximately isometric, i.e. P^Pn = I + 0{e^~^^). 

We interpret H^j-j- as the effective Hamiltonian up to order e^. The leading term 
in its symbol is hQ{y, k) = Em{k + ^{y)) + This is the well-known "Peierls 

substitution" , [20] . The symbol of hi is also quite interesting, and we discuss it in 
§4. 

Proof. As in §2 it will be convenient to work with 

H,{k) = e-'^--Hj^-- = (^-i^ - ie^ + ^ + A{y) + k^ + V{x) + W{y), 

acting on functions in B, in place of He- Note that 

H^{PN{x,y, eDy,e)u) = {27re)-^ J e'''<y-^'^/'He{k)PN{x,y,k,e)u{z)dzdk. 

The Hamiltonian He{k) can be written as He{k) = Ho{k) + eHi{k) + e'^H2{k), 
where 

Hoik) = {-id, + + k + A{y)f + V{x) + W{y) = Ho{k + A{y)) + W{y) 

Hi{k) = -2i{-id, + ^- + k + A{y)) ■ dy - idy ■ A{y), (8) 
and H2{k) = —Ay. 

We will simply construct the symbols of the pairs (/iq, -^0)5 (^15 -^1)5 • • • ) succes- 
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cancel the order zero terms in (7) we set ho{y, k) — E^aik + A{y)) + W{y) and 
Fo(a;, J/, k) = '^{x, k + A{y)). Then, since the symbol of Fo{x, y, eDy)ho{y, eDy) is 

Fo{x, y, k)ho{y, k) +j^{x, y, k) ■ ^{y, k) + 0(6% 

we must have 

(Hoik) - ho)F^ - . ^ + {H^{k) - h,)Fo = 0. (9) 

By the Fredholm Alternative in L^{E), we can solve (9) for Fi if and only if 

{Fo{;y,kl—^^{;y,k) ■ ^ + {H,{k) - h,{y,k))Fo{;y,k)) = 0, 

where (•, •) denotes the inner product in L'^{E). Hence, we choose 

h,{y, k) = (Fo(-, y, k),—-^{; y, k) ■ -^{y, k) + Hi{k)Fo{; y, k)) and (10) 

I ok ay 

F,{x,y,k) ^ {Ho{k)-hQ{y,k))-\^-^--^ + hiFo-Hi{k)Fo)+ai{y,k)Fo, (11) 

where {Ho{k) — ho{y, k))~^ denotes the inverse which maps the orthogonal comple- 
ment of Fo(-, J/, k) in L^{E) onto itself. We recall that 

Hoik) - ho{y, k) = Ho{k + A{y)) - Em{k + A{y)). 

We determine ai{y,k) by the requirement that P^Pn = / + 0{e^^^). Given 
/,5f e 5(M^), this implies 



E 



dx / [{FQ{x,y,eDy) + eFi{x,y,eDy))f][{Fo{x,y,eDy) + eFi{x,y,eDy))g]dy 



I 



fgdy + 0{e'). (12) 



Since \Fo{x, y, k)\'^dx = 1 for all (k, y), one can calculate Re{ai(y, k)} from (12), 
and sec that it is smooth and bounded, by the pseudo-differential calculus. We 

choose Im{ai(y, A;)} = 0. 

The calculation of (Fj, hj) for j > 1 proceeds in the same manner: we calculate 
the terms of order in the symbol of the left hand side of (7) which come from 
(FjT, hi) for I < j and choose (Fj, hj) so that the left hand side of (7) has the desired 
form up to terms of order e-^"*"^. At each stage we use the Fredholm Alternative in 
LF'{E), and Fj is only determined modulo a term of the form aj{y,k)Fo{x,y,k). 
Then we choose the real part of aj{y, k) so that PjPj = I + 0{e^~^^), and take the 
imaginary part of aj to be zero. Continuing in this way we complete the proof of 
the Theorem. 

Remark 2 // we set Hn = PnP^j then Hn is approximately a projection: 
= IIjv + 0{e^'^^) and IItv = H^. The equation (7) implies that 

n* TJ TJ r>* I /nr ^N+l\ 
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Hence 

If we replace P^Hf^ff in the equality above by the left hand side of (7),we obtain: 

PnPnHs = HsPmPn + C(e^+i). 

Thus Hn is a projection which commutes with to order 0(e^^^) as in [19]. 

In [19], the construction of the almost invariant subspaces is based on the method 
of Nenciu-Sordoni [17] and Sordoni[18] (see also [10], [14], [15]). This method is 
heavily related to the construction of Moyal projections. 

4. Relations with Previous Work 

To relate the resuhs here to what has aheady been done we need to complete 
the calculation of the effective Hamiltonian up to terms of order e^, i.e. to compute 
the symbol hi{y,k) from (10). The explicit computation of 

(Fo(-, y, /c), ^i(/c)Fo(-, y, k)) ^def Ei{y, k) 

is contained in the computations in [5] (it is also in [8] with a small error - see 
Remark 1 in [5]). When one replaces k{y,s) and \l/(x,y,s) in [5, pp. 7601-3] by 
k-\-A{ii) and ^'(a;, k-\- A{ii)) respectively, Ei^y., k) is "i/i" in the notation of [5, (25)] 
and we have 

1 d dE^ik) ~ d^j^. dEmjk) 

where k = k + A{y). Here B{y) = V x A{y) and L = 

L = Im ((Mto, *) g^), (Mfa, k)—. — ),(Mfa, k)—. —) 

with M{y,k) = IIo{k) — ho{y,k). The vector L is an angular momentum and 
L ■ B contributes the "Rammal- Wilkinson" term to the energy, cf. [1] . Adding the 
additional term from (10) to Ei to obtain /ii, we obtain 

^^ = ki§y- ^^"^al^^^^^ - ^ • ^ - ^ + ^(^)' ^ + ^(^)))' (13) 

where 

and y and k are defined by the Hamiltonian system 

. ^ d{E,r.{k + A{y)) + W{y)) d{E^{k + A{y)) + W{y)) 

^ dk dy ' 

Thus one recognizes i(\I/(-, k+A{y)), ^'(•, k+A{y))) as the term generating the Berry 
phase precession, cf. [13], [22]. Comparing (13) with [19, (22)] (in the case / = 1), 
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the choice of sign in the magnetic potential, A{y), and the use of Weyl quantization 
in [19]. The sign of the Berry phase term in (13) may appear inconsistent with 
[5, (29)], but it is not. In [5] (y. A;) was the vector field from the Hamiltonian 

-Emik + A{y)). 

In [5] and [8] instead of introducing effective Hamiltonians we constructed wave 
packets. These packets are nonetheless related to effective Hamiltonians in that 
one can compute what the effective Hamiltonian must be - assuming that there is 
one - from the packets. To see this one can proceed as follows. The packets have 
the form (here s = et and W{y) = 0) 

u{x, y, s, e) = e^^(^'^)/^[/(2/, s)^{x, ^ + A{y)) + 0{e)], 

where </> and / are solutions of 

Here all fimctions of (/c,y) are evaluated at /c = /c(y, s) = dyCp, and D = {l/2)dy ■ 
{dkEm{dy(j) + A{y))). Assuming that the evolution of / is governed by an effective 
Hamiltonian i^e// = ho{y, ^Dy) + ehi{y, eDy) + C(e^), we must have (on bounded 
intervals in s) 

[gi.He/,Ag^0(.,o)A^(.^ 0)](2/, s) = e'^^y''^/'f{y, s) + ^(e^). (15) 
Differentiating (15) with respect to s, one concludes 

'-H,ff{e'^(y'^yj{y, s)) -C-^f + |^)e^^(^'^)/^ + 0{e). (16) 
Using the symbol expansion from the pseudo-differential calculus 

e->1'''H,,f{e'^-''f) = ho(y,k)+ 

^' ■ I w 4 Sk<^- *>5S;'^' ^'i^(^) (1^) 

Substituting (17) into (16) and comparing the result with (14), one recovers the 
formulas given earlier for /?.o(y, k) and hi{y, k). 

In [8] we were unable to reconcile our results with those of Chang and Niu, see [3] , 
[4] and also [23]). We thought that this might have resulted from different choices 
of scales. This is partially true, since Chang and Niu do not distinguish the scale 
y = ex, but, as Panati, Spohn and Teufel point out in [19], the differences largely 
disappear when one considers the Heisenberg formulation of quantum dynamics. 
Letting a{y, eDy) be an observable, the propagation of a in the Heisenberg picture 
is given by 

a{s) = e-''^^ff/'ae''"-ff/' or = i[a, H^ff]. 

If one considers this propagation at the symbol level, then the symbol of a is 
propagating along the trajectories of the Hamiltonian system 
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Hence one can consider the contribution of hi as an order e correction to the classical 
Peierls dynamics arising from ho. This is the point of view taken in [19]. However, 
it is worth noting that the wave packets are propagating along the trajectories from 
ho with a precession in their phases arising from (the imaginary part of) hi. 
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